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Abstract 

The old problem of a singular, inverse square potential in nonrelativistic quantum mechanics 
is treated employing a field-theoretic, functional renormalization method. An emergent contact 
coupling flows to a fixed point or develops a limit cycle depending on the discriminant of its 
quadratic beta function. We analyze the fixed points in both conformal and non-conformal phases 
and perform a natural extension of the renormalization group analysis to complex values of the 
contact coupling. Physical interpretation and motivation for this extension is the presence of 
an inelastic scattering channel in two-body collisions. We present a geometric description of the 
complex generalization by considering renormalization group flows on the Riemann sphere. Finally, 
using bosonization, we find an analytical solution of the extended renormalization group fiow 
equations, constituting the main result of our work. 
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I. INTRODUCTION 



Exactly solvable problems play an important role in physics. They provide a backbone of 
our understanding and allow to develop perturbation theory on the basis of an exact solu- 
tion. In nonrelativistic quantum mechanics the harmonic oscillator and the Coulomb central 
potential are certainly the two most important examples. Another prominent example of 
an exactly solvable problem in quantum mechanics is the central, inverse square potential 

V{r) = (1) 

Most remarkably, in any spatial dimension the potential ([1]) is classically conformal invariant 
because it is a homogeneous function of degree -2 and has the same scaling as the nonrel- 
ativistic kinetic energy. Hence, the classical action is invariant under the nonrelativistic 
scaling r — pr, t — > p^t. The quantum problem with the potential (II]) has a lone history 
and is discussed both in standard textbooks 1, 2| and in the scientific literature jsl, 0]. More 
recently there were a number of renormalization group (RG) studies of this problem 



15|. 



The quantum physics of the inverse square potential is well understood. The interaction 



dl]) is an example of a singular potential 16| and must be treated with care. It is known that 



18|. The anomaly has its origin in 



for the repulsive and weakly attractive coupling (k < HcrY, the scale symmetry is preserved 
at the quantum level and the theory provides an example of nonrelativistic conformal field 
theory [17|. On the other hand, for strong attractions {k, > Kcr) a discrete, geometric bound 
state spectrum develops in the two-body problem, and the continuous scale symmetry is 
broken to a discrete subgroup by a quantum anomaly 
the singular short- distance behavior of the inverse square potential. In the RG language the 
anomaly manifests itself as a limit cycle instead of a scale-invariant fixed point. 

The inverse square potential is a paradigmatic system for nonrelativistic conformal in- 
variance and scale anomaly. Remarkably, there is a number of different physical systems, 
which are described (often only in some restricted domain) by the inverse square potential: 



The celebrated Efimov effect, first derived in [19|, consists in the formation of a tower 
of three-body bound states of identical non-relativistic bosons^ interacting through a 
short range potential. Exactly at resonance (unitarity regime) all scales drop out of the 



^ Kcr represents the critical attractive coupling. In d spatial dimensions it is given by Kcr = ^'^ ■ 

^ More generally, for the occurance of the Efimov effect it is sufficient that at least two two-body subsystems 
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problem and the three-body spectrum becomes infinite and geometric. The three-body 
)roblem in quantum mechanics is treated most easily in hyperspherical coordinates 



20|. Employing the adiabatic hyperspherical approximation 2l|, Faddeev decompo 



sit ion of the wave function 



22| and restricting to the zero total angular momentum 



sector, one arrives at the remarkably simple one-dimensional effective equation 23] : 



1/4' 



ip{r) = Eip{r) 



(2) 



with the Efimov parameter sq ~ 1.0062. This is a one-dimensional radial Schrodinger 
equation with the inverse square potential in the overcritical regime (k > j). 

The interaction of a polar molecule with an electron in three spatial dimensions can 
be approximated by a point dip ole- charged particle anisotropic potential 

cos 6* 



V{r) 



(3) 



where the angle 6 is measured with respect to the direction of the dipole moment. 



As was demonstrated in 



2J], this anisotropic conformal potential can be reduced to 



the effective isotropic inverse potential ([T]) in the zero angular-momentum channel. 
Remarkably, the dipole-electron system exhibits all interesting features characteristic 
for the potential ([I]). 

At the critical coupling the theory undergoes a transition from the conformal 
(k < Kcr) to the non-conformal (k > Hcr) regime. It has been demonstrated recently 
in 15| that this transition is closely related to the classical Berezinskii-Kosterlitz- 
Thouless (BKT) phase transition in two dimensions. In particular, the energy B of 
the lowest bound state near the critical coupling in the non-conformal regime 
vanishes like 

5 ~ exp ( ^ ^ ) (4) 



which is analogous to the behavior of the inverse correlation length as the BKT tran- 
sition temperature is approached from above 25||. 



have s-wavc bound states close to the zero-energy threshold unless two of the three particles are identical 
fermions. 
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Quite unexpectedly, the inverse square potential arises in the near-horizon physics 
of some black holes. More specifically, this is the case for a massive, scalar field 
minimally coupled to gravity on the non-extremal spherically-symmetric Reissner- 
Nordstrom (RN) black hole background. It was demonstrated in 



that the Klein- 



26| 



Gordon field equation of the scalar field reduces to the effective Schrodinger equation 
with the overcritical (k > inverse square potential in the near-horizon limit. The 
related problem of a scalar particle near an extremal RN black hole was treated in 
leading to a similar finding. 

Finally, the conformal inverse square potential appears naturally in the context of 
vacuum AdS / CFT correspondence . The field equation for a scalar field of mass 
m in the Euclidean AdSd+i is 

- ^dr<P - ^0 - = 0, q' = {qy + (f, (5) 

where we transformed to the momentum space on the AdSd+i boundary [x^, x) — >• 
q), and r denotes the radial direction in the AdSd+i space. We can change variable 
(p = r^^~'^y^^ and obtain 

which is a one-dimensional Schrodinger equation with an inverse square potential 



15|, 



of strength k = — — and energy E = —q^. As was emphasized in 
the overcritical coupling k > n^r corresponds to the violation of the Breitenlohner- 
Freedman bound in AdSd+i- 

In this work we study the quantum problem of the inverse square potential using a functional 
renormalization method. The plan of this paper is the following: In Sec. |TT]we introduce our 
method, the physical system of interest, and derive the fiow equations using the sharp cut-off 
regulator. A general mathematical discussion of the flow equation is performed in Sec. IIIII 
Our central result is derived in Sec. IIV[ where we extend the analysis to the complex plane, 
discuss the flxed point structure and flnd a numerical solution of the extended set of flow 
equations. Additionally, we provide a physical interpretation, geometric description and 
motivation for the complex extension. In Sec. |V]we bosonize the interaction. This allows us 
to view the problem from a different angle and, most remarkably, obtain an analytic solution 



of the generalized complex flow equations. We draw our conclusions in Sec. IVII Finally, in 
Appendix |X] we present the renormalization group flows on the Riemann sphere. 



II. THE METHOD, THE MODEL AND THE FLOW EQUATION 

In this work we use the functional renormalization group and calculate a scale-dependent 



effective action functional 28| (for reviews see 29l-l3l[|) called also average action or flow- 
ing action. The method is formulated in Euclidean spacetime and employs the Matsubara 
formalism if one works at finite temperature. The flowing action F^ includes all fluctuations 
with momenta q ^ k. In the infrared (IR) limit k ^ the full quantum effective action 
F = Ffc^o is recovered and the problem is solved. In practice the dependence on the scale 



k is introduced by adding a regulator term Rk to t 
action F^ obeys the exact functional flow equation 



le inverse propagator F(2). The flowing 
281 ]. and its bosonic version is given by 



dkTk = ^TrdkRk (rf^ + Rtr' = ^Trdk In(Ff) + R^). (7) 

The functional differential equation for F^ must be supplemented by an initial condition 
Ffc^A = S. The "classical action" S describes the physics at the microscopic ultraviolet 
(UV) scale k = A and is assumed to be known. In Eq. (I7j) Tr denotes a trace which sums 
over momenta, Matsubara frequencies, internal indices, and fields. The second functional 

(2) 

derivative F^ is the full inverse field propagator, which is modified by the presence of the 
IR regulator R^. The momentum dependent regulator function Rk{q) must satisfy three 



important conditions, but otherwise can be chosen arbitrary |28|. The regulator dependence 
of the flowing action drops out for A; — 0, and one obtains the exact solution of Eq. ([7]). In 
the second form of the flow equation ([7]) dk denotes a scale derivative, which acts only on the 
IR regulator Rk- This form is convenient because it can be formulated in terms of one- loop 
Feynman diagrams. It is also convenient to introduce the RG "time" t = ln{k/ko), where ko 
is an arbitrary reference scale. The RG evolution towards the UV (IR) corresponds to t > 
(t < 0). In the following we will use both t and k. 

In most cases of interest Eq. ([71) can be solved only approximately. Usually one adopts 
some type of expansion of F^ and then truncates at finite order. This leads to a finite 
system of ordinary differential equations. The expansions do not necessarily involve any 
small parameter and are generally of non-perturbative nature. 
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Although our method allows us to address the quantum many-body problem at finite 
temperature, in this work we are interested only in the few-body (vacuum) physics, which 
is characterized by vanishing density (n = 0) and vanishing temperature (T = 0). In this 
case the effective action Tk=o, being the generating functional of the IPl vertices, can be 
easily related to the different scattering amplitudes and bound state energies. Additionally, 
numerous simplifications arise in the structure of the fiow equations (for a detailed discussion 



see 



32 



Finally, we note that the three-body quantum problem for fermions and bosons was 



treated with functional renormalization in 
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35|. In this work we apply functional renor- 



malization to solve an arguably simpler problem, and our aim is twofold. First, the inverse 
square potential is a paradigm and it is important for understanding of other, more challeng- 
ing problems. Second, working with this simple system, we develop a new renormalization 
group method of complexified fiows in this paper. The method may be useful for more 



technical few-body problems in the future 36|. 



In this work we study the nonrelativistic quantum mechanical problem of identical bosons 
interacting through a long-range potential ([T]) in (i spatial dimensions. The many-body field 
theory is defined in the UV by the microscopic action Se 

j-l/T I- 

Se[(P, iJ*] = dr d'^xip^T, x)[dr - A - /i]V^(r, x) 
A /"^/^ 

/ dTd'^X^P*{T,x)tP*{T,x)llj{T,x)ij{T,x) (8) 



2 

dr / d'^xd%ij*{T,x)i/j*{T,y)- —^lj{T,y)tlj{T,x). 

J 



Our convention is h = 2M^ = 1 with the boson mass M^. We work in the Matsubara 
formalism with Euclidean time r G (0, 1 /T). In what follows we will be interested exclusively 
in the few-body (vacuum) physics. The vacuum state is characterized by zero density, which 
corresponds to zero chemical potential (/i = 0), and zero temperature (T = 0). The bare 
action ([8]) is invariant under a global U{1) transformation and possesses Galilean space-time 
symmetry. The microscopic bare parameter k characterizes the strength of the long-range 
potential and is positive in the attractive case. We augmented the theory by a four-boson 
contact interaction term with a coupling Aw,. This is a consequence of the fact that the 

nn 

inverse square potential is singular at the origin p|, |16[. On its own the singular potential 
is not sufficient to define a quantum mechanical problem and must be augmented by the 



boundary condition at the origin [l5j. We will see that the introduction of in the UV 
determines the r = boundary condition and makes the problem well-defined. 

For our method it is convenient to switch to momentum space. The Fourier transform of 
the potential in d spatial dimensions reads 



(4vr)'^/2r(c//2 - 1)1/1^-^ 2<d<5. 



(9) 
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The restriction to dimensions d in the range 2 < d < 5 can be understood easily by the 
fact that on the one side, d = 2 is a natural lower dimension, in which the integral is IR 
logarithmically divergent. On the other side, the upper bound d = 5 can be relaxed, if we 
modify the Fourier integral by the introduction of a UV suppression factor exp(— e|r|) and 
perform the limit e — )■ in the very end. Hence, we will use 



{47iY/^T{d/2-l)\l 



2~d 



FS) = ^ — d>2. (10) 

In what follows we consider only d > 2. First, we consider a momentum-independent 
(pointlike) truncation for the fiowing action 

tp -t- r^^ 

Q 



K j FSW{Qi)i^{Q2W{Q^)'^{QA)K-Qi + Q2 - Qs + Q4) (11) 

J Qi .0.7. Qa 



T 



r{QM{Q2W{Q^)^{Qi)K-Qi + Q2 - Q3 + Qa] 

Qi,Q2,---,Q4 



where Q = {uj,q) and = J ^ J The vector I = - qi = - q^^ gives the 

spatial momentum transfer during a collision and / = |/|. In the nonrelativistic vacuum the 
propagator of the elementary field ip is not renormalized because the only diagram, which 
contributes to its fiow, contains a hole (antiparticle) in the loop. As there are only particles 
but not holes in the nonrelativistic vacuum the propagator keeps its microscopic form. 
Remarkably, the coupling k, characterizing the strength of the long-range potential, is 
also constant during the renormalization group fiow (we discuss this issue in more detail in 
Sec. |V]). The only coupling which fiows during the RG evolution in our truncation is the 
contact coupling A^. Its fiow equation is 

* ^ A \iuj + l^ + Rk{L)){-iu + P + Rk{-L)y ^ ' 



Figure 1: Flow equation of the four-particle contact coupling in form of Feynman diagrams. The 
solid lines with an arrow denote the boson tp regularized propagator, the dashed lines correspond 
to the long-range interaction vertex, and the dark blobs represent the contact coupling A^. 

The flow equation is depicted in terms of Feynman diagrams in Fig. [1] In order to proceed 
further we must specify the cut-off function Rk{L). 

Nevertheless, we note that already at this point it is possible to identify the generic form 
of the flow equation which reads 



dtX^ = + PX^ + 7 



(13) 



with a, /3, 7 G M, which depend on the coupling k, dimension d and a concrete choice of the 
cut-off function. We postpone the general analysis of equation (fT3l) until Sec. IIIII 
It is convenient to rewrite Eq. (|T2l) as 

(14) 



dtX^ = JqaXI, + 4 Ji_dKA^, + 4 J2,dK^ 



where we defined the cut-off dependent integrals 

duj dH 



-A 



(15) 



2n {2nY * {luj + P + Rk{L)){-iuj + P + Rk{-L)) 

where F2{1) = Fd{l)"'. In the rest of this section we finish the computation by specifying 
the sharp regulator Rk{L) = {iu + /^) {^jfjTzj^ — ij- 

The cut-off function Rk{L) = {icu + /^) ^g(|2zpy — cuts off quantum fluctuations 
sharply, i.e. it totally suppresses the modes with < /c^ during the renormalization g roup 



evolution. This type of cut-off was used in the closely related three-body problem in 
The cut-off satisfies 

9{l^ - 



33|. 



iio + P + Rk{L) iu + P' 

Employing this property we can calculate Jn,d defined in Eq. ( ITSll explicitly as 

du dH 



(16) 



J, 



n,d 



27r (27r) 



TrSd 



m) 

[lu + P){-iuj + P) (27r) 



-,k''-'F2{k) (17) 
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with Sd = r{d/2) being the area of a unit sphere in d- dimensional space. The resulting Jn^d 
can be readily substituted in Eq. ( fT4l) 

"^'^"^ - 'jA^^wwrn " ^ " ~ (rf-2)2 ■ ^^^^ 

The contact coupling has a naive (Gaussian) scaling dimension [A^] = 2 — d and is IR 
irrelevant in > 2. Our aim is to investigate the fixed point structure of Eq. f|T8l) . For 

V« — (47r)'*/2r(d/2) 



this reason we introduce a dimensionless, rescaled coupling A^r = , \d/2ruii'}\ ^^"^ 



equation reads 

5*A^« = -AJ^ + (^-^ + - 2^ A^« - j^^r (19) 

As will be shown in Sec. IIIIl the discriminant D of the quadratic /3-function determines the 
overall behavior of the solution. In our case D = —Ak + (d — 2)^. The critical Kcr is defined 
by the condition D = yielding 

^cr = (20) 

which is in agreement with the quantum-mechanical non-perturbative calculation 15[. For 
K < Kcr (weak attraction and repulsion) the dimensionless coupling A^^ has two real fixed 
points 

x± d-2 , Ud-2f 

2^—^ r-T^-^ 

and the theory is scale invariant in the IR respectively UV. For k > (strong attraction) 
the coupling ceases to have real fixed points and scale invariance is lost. 



III. GENERAL ANALYSIS OF THE FLOW EQUATION 

The fiow equation for found in Sec. Hi] has the general form 

^A^fi(t) = aX^nitf + p\^R{t) + 7, (22) 

where Q!,/3,7 G M are numerical coefficients. Without loss of generality we consider a < 0. 
This choice corresponds to the result obtained in Sec. [Tll The form of the solution of Eq. 
(122|) is determined by the sign of the discriminant D = fS"^ — 4a'-f of the quadratic /3-function. 
We consider the three different cases depicted in Fig. [2] 

• D>0 
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Figure 2: The /3-function of the coupHng A^r for a = —1 and /3 = 0. The three hnes correspond 
to different values of 7 with D = 1 (solid blue), D = (dashed red) and D = —1 (dashed black). 

In this case the /3-function has two real fixed points and A^^ with A^^ < A^^. 
The solution of Eq. ( 122|) depends on the interval, where the UV initial condition A^j^ 



belongs. For the initial condition A:^^ 

-(3 - y/D tanh 



< A^i? < A^^ the solution is attracted towards 



the IR by the fixed point A^,^ and has the form: 



2a 



(23) 



with rj fixed by the UV initial condition X^R{t = 0). For the initial condition X^r > 
Xuv the flow runs into a positive divergence (Landau pole), but reemerges at negative 
infinity and subsequently approaches A^,^ in the IR. The explicit solution in this case 



reads 



-/3 + coth 



i-t + v) 



2a 



(24) 



with ?7 fixed by the UV initial condition. In fact, the case k = 0, which leads to D > 0, 
corresponds to the well-studied case of a contact interaction in atomic physics. In this 
context the appearance of the Landau pole signals the presence of a weakly bound 
molecular (dimer) state. Finally, for the initial condition X^r < A^^ the RG flow is a 
smooth, monotonic function X^R{t), which approaches A^^ with the explicit solution 
given agam by Eq. ([M]). Notably, for D > the IR value A^r(A; = 0) = A^^ is not 
sensitive to the concrete choice of the initial condition. We offer an elegant geometric 
description of the flow in Sec. |IV] and Appendix 1X1 
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Figure 3: The renormalization group flows of the couphng A^r for a = —1 and /3 = and initial 
condition A^/j(i = 0) = 0. The different lines correspond to three different choices of 7 with 
D = l,r] = (solid blue), D = 0, 77 — )■ 00 (dashed red), and = — 1, = (dashed black). 

This is a limit of the previous case {D > 0) when the two fixed points merge. Excluding 
the trivial choices a = /3 = 7 = and a = /3 = 0;77^0, the /3- function has a single 
fixed point A^^ — (^^^ Fig. The RG equation takes the form 



dt 



with the solution 



A^ij(t) — A!,^ - 



at + Tj 



(25) 



(26) 



where 77 is fixed by the initial condition. The running of the coupling \^,R{k) is log- 
arithmic (which corresponds to the marginal deformation) and for 77 < it hits a IR 
Landau pole. The coupling runs into the pole at the scale t = nevertheless the 
RG evolution can be extended beyond this scale and approaches X*^p^ at /c = 0. 

D < 

In this case there are no real fixed points. The formal solution can be written as: 



Xil>R(t) 



-P + v^^tan 



-D 



(t + v) 



2a 



(27) 



where rj is fixed by the initial condition. This solution is periodic with a period 



2tt 



-D' 



Remarkably, in this case the UV initial condition \^R{k = A) plays an 



important role as it determines the infrared value \^R{k = 0)." 



^ This is analogous to the necessity to introduce the well-know three-body parameter in the context of the 
Efimov effect [23 1. 
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We plot the RG flows of X^n in the three cases in Fig. [31 While for D > the flows 
approach a fixed point in IR, the renormalization group evolution develops a limit cycle^ 
for D < 0. As the RG scale k can be related to the energy of the particles, the physical 
interpretation of the limit cycle solution is clear: During the RG flow one hits bound states, 
manifesting themselves as divergences of the coupling X^r. And since there are inflnitely 
many divergences one has an infinite tower of bound states with a geometric spectrum for 
the case D < 0. 



IV. COMPLEX EXTENSION 

In this section we discuss the different subcases, considered in Sec. IIIH in more detail 
and extend the analysis to complex values of the interaction coupling. We also present a 
physical interpretation of this extension. 



A. Negative discriminant: complex fixed points 

The renormalization group flow equation in our point-like approximation is^ 

dtX = aA^ + /5A + 7 (28) 

with a, /3, 7 e M, and one can vary the parameters /3 and 7 by considering different n and d 
in Eq. ( IT9l) . It is sometimes useful to express Eq. (l28l) in the alternative form 



a^A = a(A - A,)2 + A7, ^* = -^' = ^^^^ 

If the discriminant D = P'^ — 4a7 = —4aA^ < 0, the /3-function has a pair of complex roots 

A^^-^^'^ = A.±.^g3. (30) 

2a a 



^ Strictly speaking, to find a limit cycle one needs at least two couplings connected by the RG flow equations 
(see, for example, Eqs. ((63l l65|) ). In the case of a single coupling constant an infinite (unbounded) limit 
cycle appears only if there are periodic real discontinuities in the RG flow. We provide an elegant 
description of the infinite limit cycle on the Riemann sphere in Sec. IIVI 

^ This equation coincides with Eq. ([22]) ■ To simplify notation, we denote the coupling A.^/? as A in this 
section. 
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For this reason it is natural to consider the RG evolution of a generally complex variable 
A = Ai + 1X2 in Eq. fl28l) . The resulting flow equations in the complex plane now read 

dtXi = aXl — aXl + /3Xi + 7 

(31) 

dtX2 = 2aAiA2 + 13X2. 
Notably, if we start on the real axis (i.e. X2{t = 0) = 0), the flow remains real 

dtXi = aXj + /3Ai + 7 

(32) 

dtX2 = 



with the periodic solution 



2a 



where rj is flxed by the initial condition Xi(t = 0). The solution can be nicely identified 
with a limit cycle if we map the complex plane of the Riemann sphere. On the Riemann 
sphere the real line corresponds to a great circle and the solution ( l33l) traverses this circle 
periodically (see Fig. S] (A) and Appendix Remarkably, on the complex plane the real 
axis is the separatrix of the two complex fixed points. 

Let us investigate the properties of the fixed points A^ by considering the stability ma- 
trices 

MS ^ ^U.. (34) 

A straightforward computation shows 

^ ( -b^\ ^ , 

M± = , = ±i^/\D\. (35) 

The stability matrices have a pair of complex conjugate pure imaginary eigenvalues 

nf = ±iyW\ 4 = {nfy. (36) 

Therefore, the local fiow near the fixed points has a form of a circle. The sign of 6^ 
in the matrix (153]) determines the orientation of the circle. For b^ > the circulation is 
anticlockwise, while for 6^ < it is clockwise. Hence the fixed points A^ have the opposite 
circulation. 

The phase portrait of the RG fiow in the complex plane, computed numerically for specific 
values of parameters, and the example of the RG fiow are depicted in Fig. [51 We observe 
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A = oo A = 00 




Figure 4: The real flows of A projected onto the Riemann sphere: Great circle with arrows pointing 
towards the UV represents the real line. (A) For negative discriminant the fixed points A"*" and A~ 
are complex, and the flow traverses the real great circle periodically generating a limit cycle. (B) 
For positive discriminant the IR fixed point A^ and the UV fixed point A~ lie on the real great 
circle. The UV fixed point can be reached from the IR fixed point via two different paths- the 
"short" path I and the "long" path II. 



Ai 




Figure 5: (A) The phase portrait of the flow equations (j3ip for the specific choice of parameters 
a = 7 = —1 and /3 = 1. Arrows denote the direction towards the UV. (B) Corresponding periodic 
flows of the real part Ai (blue) and the imaginary part A2 (red) of the complex coupling A. 

that the periodic divergences of the pure real solution (l33l) are regularized and the flow of 
the imaginary part A2 develops a tower of resonances. The analytic solution of Eq. fl3T]) for 
-D < is given in Sec. |Vl 
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B. Positive discriminant: real fixed points 



If the discriminant D is positive, the /5-function has a pair of real roots 



-/3±\/D 
2a 



(37) 



As before, we generahze the couphng to complex values A = Ai + 1X2 and obtain the pair 
of coupled differential equations (13T|) . The phase portrait and a specific solution of the flow 
equations fl3T]) for D > are depicted in Fig. [61 We find the analytic solution of Eq. f l3T|) 
for D > in Sec. |Vl 

The character of the fixed points A-t can be determined from the stability matrices M^, 
defined by Eq. (IMl) . For D > 0, we obtain 



The sign of the real eigenvalue determines whether the fixed point is UV attractive or 
repulsive. The eigenvalue is positive and hence the fixed point A+ (left fixed point in 
Fig. Et^A)) is UV repulsive, meaning that as the sliding scale k is increased the flow is driven 
away from A''". On the other hand k,~ is negative and therefore the fixed point A~ (right 
fixed point in Fig. [6](A)) is UV attractive. Notice that any two-component vector is an 
eigenvector of the fixed points A^. 

We propose a geometric interpretation for the different behavior of the real solution 
in dependence on the initial conditions, which we observed in Sec. IIIII For a positive 
discriminant D both the UV fixed point A~ and the IR fixed point A"*" are situated on the 
large real circle on the Riemann sphere (see Fig. H] (B) and Appendix We note that 
the UV fixed point can be reached from the IR fixed point following two different paths. 
Depending on the initial conditions Ain = A'^ ± e, the flow can follow either a "short" path 
I or a "long" path II (see Fig. H] (B)). Path I corresponds to the solution (l23l) . which is a 
standard way how to regularize the inverse square potential at k < On the other hand, 
the path II traverses A = 00 on the Riemann sphere and corresponds to the solution flM]) 




(38) 



with the degenerate eigenvalues k- 




(39) 
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Figure 6: (A) The phase portrait of the flow equations (I31|) for the specific choice of parameters 
a = — 1 and /3 = 7 = 1. Arrows denote the direction towards the UV. (B) Corresponding 
renormahzation group flows of the real part Ai (blue) and the imaginary part A2 (red) of the 
complex coupling A. 



with the Landau pole. The divergence during the RG evolution is identified with a single 



bound state, which, in the case o 
extensively in atomic physics 



K = 0, becomes the the well-known shallow dimer studied 
23| . We notice that on the full Riemann sphere the path II 
can be continuously deformed into the path I. Thus, the introduction of a small imaginary 
initial condition \2{t = 0) for the complex extended flow, which leads to a small deformation 
of the path II, regularizes the divergence in the real part of the coupling (see Fig. [6] and 
Appendix lAl). 



C. Zero discriminant: degenerate fixed points 

Finally, for Z) = 0, the pair of roots of the /3- function becomes degenerate 

A, = A+ = A- = (40) 

which corresponds to a single, real- valued fixed point A*. 

As it turns out, the stability matrix Mij = ^\\, vanishes in the case of D = 0. This 
corresponds to a logarithmic (marginal) renormahzation group flow in the vicinity of the 
fixed point and local properties of the degenerate fixed point A* are accounted for by the 
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Ai 




(A) (B) 

Figure 7: (A) The phase portrait of the flow equations (I31|) for the specific choice of parameters 
a = —1, P = 2 and 7 = — 1. Arrows denote the direction towards the UV. (B) Corresponding 
renormalization group flows of the real part Ai (blue) and the imaginary part A2 (red) of the 
complex coupling A. 

second derivative of the /3-function Kiji at A^, 

The only nontrivial components are 

Kiu = -Ku2 = K212 = K221 = 2a. (42) 

Again, we computed the phase portrait and a specific solution of the flow equations fl3T|) 
and show the results in Fig. [71 We observe that the Landau pole of the real part of A is 
regularized and the imaginary part of A develops a single resonance. 

D. Physical interpretation of the complex extension 

One may wonder about the physical meaning of the complex contact coupling A intro- 
duced in this section. To answer this question, we consider inelastic scattering^ of particles 

^ By definition inelastic collisions change the internal state of the colliding particles and, henee, at least one 
of the colliding particles must be composite and have some internal structure. For example, in experiments 
two atoms might collide and fall into energetically deeper lying internal states. The excess in energy will 
be converted into kinetic energy and the atoms will be lost if the released energy suffices to overcome the 
trapping potential. 
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in a central potential in quantum mechanics. Following [2[ in three spatial dimensions the 
scattering amplitude in the center of mass frame can be expanded in partial waves 

Z=oo 

Me) = 5^(2/ + i)Mk)Pi{cos{e)), (43) 

where Pi{cos6) are the Legendre polynomials. From the unitarity constraints the partial 
waves scattering amplitudes fi{k) as functions of scattering momentum k are given by 

- (^^) 

Here, gi{k'^) is an even function of k which is generally complex in the case of presence of 
inelastic channels. 

For low-energy scattering of particles, interacting through a short-range potential, only 
the s-wave contribution / = /o is substantial and gik"^) = goik"^) can be expanded as 

g{e) = -a~' + ^r,sk' + ..., (45) 

where in the case of elastic scattering a and are real and denote the scattering length 
and effective range. However, in the inelastic case a may also be regarded as a complex 
scattering length a = a + z/3 and the effective range r^s might also be complex js^j. In the 
forthcoming argument it is sufficient to consider low energy scattering and therefore we keep 
only the scattering length a and neglect the second term in Eq. fHSjl . 

The generalized optical theorem for the scattering of indistinguishable particles 

k k 

lmf{k) = -^(Ttot = -^((^el + CTin) (46) 



holds also in the case of general inelastic scattering |2[. In this case the total scattering 
cross section atot is the sum of the positive elastic aei and inelastic crj„ contributions. The 
positivity of (Tj.„ and the optical theorem P6|) imply that the imaginary part of the scattering 
length P must be negative. At low energies and for sufficiently short-range interactions the 
complex scattering length fully determines the elastic and inelastic cross sections, and for 
indistinguishable particles the result reads ^ 

Gel = Sn\a\^{l-2k\(3\) 
a,„ = 87r^(l-2A;|/3|). 
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We observe that a non-vanishing imaginary part of the scattering length (3 is required to 
obtain an inelastic contribution to the total cross section'', which results in loss of particles. 

The presented arguments can be connected to a complex generalization of A by the simple 
observation that in a system of identical bosons with short-range interactions the physical 
scattering length a is related to the IR value of the dimensionful, contact coupling in 
d = 3 via the simple formula (in the case of k = 0) 

X^{k = 0) = -Svra. (48) 

Thus, allowing for complex values of A^ during the RG flow corresponds to the presence of 
inelastic two-body collisions leading to particle loss. This provides a physical interpretation 
of the complex coupling A^. We must stress, however, that the local {ip'^ipY operator with a 
negative imaginary coefficient can describe inelastic scattering only to deep energetic states, 
i.e. the gap energy Egap of the state must be large compared to other energy scales in the 



problem [381] • On the other hand, decay into shallow bound states can be described by 
introduction of a non-local operator with a complex coefficient. In fact, a large imaginary 
part of the two-body interaction might also be useful in the context of atomic physics. This 
case has recently been studied for cold atoms confined to motion in ID. In particular it 
has been shown that a Tonks-Girardeau gas, where ID strongly repulsive bosons exhibit 



fermionic like behavior, can also be induced by the effect of strong dissipation 39|, |40| . 
meaning large imaginary two-body interaction A^2- 

Finally, we note that our extension of the RG equation to the complex plane can be used 
as an efficient numerical tool. As was demonstrated in Sec. Illlt the flow equation fl22l) has 
the periodic solution ( !27l) in the non-conformal phase. The running coupling A^,r diverges 
periodically during the RG evolution, making the numerical solution impossible beyond the 
flrst divergence. We showed in this section that adding a small imaginary part to the initial 
condition A^/j(t = 0) makes the solution of Eq. (!22l) numerically feasible. Physically, this 
corresponds to converting stable bound states to long-lived resonances. As shown in the 
next section, the problem considered in this work can be treated analytically and hence the 
numerical treatment is not necessary. Nevertheless, the extension turns out to be useful in 



^ We note that according to Q Eq. (|47p is valid for sufficiently fast decrease of the interaction potential 
at large distances (at least Thus, it is not strictly applicable in our case of the long-range 

potential, and the presented argument has to be taken as heuristic. 
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some other few-body problems, which can not be solved analytically 



36| 



V. BOSONIZATION 



In this section we approximate the four-particle vertex by an exchange of a composite 
particle. In the literature 3l| this is known as bosonization procedure and it is a powerful, 
physically transparent concept. In this work we will bosonize in two distinct ways. The 
first bosonization corresponds to the exchange of a massless non-dynamical particle in the 
t-channel. In this realization the composite particle mediates the long-range interaction 
combined with the contact four-particle interaction. Alternatively, one can approximate the 
contact vertex by the exchange of a massive particle in the s-channel. While at low energies 
and momenta bosonization simply reproduces the pointlike (momentum-independent) ap- 
proximation of Sec. ini it resolves some of momentum structure of the interaction at higher 
energies. 



A. t-channel bosonization 



In this subsection we approximate the total four-particle interaction vertex by the ex- 
change of a composite massless particle x i^i the t-channel. 




Galilean symmetry of the nonrelativistic vacuum implies that the inverse propagator P^{u,q) 
of the massless particle depends only on the absolute value of the spatial momentum q, but 
not on the frequency u. Thus, the composite x is ^ot a dynamical particle and serves only 
to mediate the two-body long-range interaction. 
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Mathematically, bosonization is implemented by means of the Hubbard-Stratonovich 
transformation. At T = /i = one adds an auxiliary part to the microscopic action ([8]) 



SE,au.[i^.r.x] = ^ [ Pxil){xi-L) + -^r{QMQ2)}y< 

^ JQuQ2,-,Q4 2P^(0 

{x(L) + -^riQ3)HQ^)}x (49) 

S{-Qi + Q2-Q3 + Q4), 

where L = {ui, I) = Q2 — Q1 = Qz — Qi is the momentum transfer (i+l-vector and x is a real- 
valued field. The modified action S'^ = SE+SE,aux defines an equivalent quantum mechanical 
problem as the action Se because the functional integral in the field x is Gaussian. 

Notably, at the UV scale we can achieve the cancellation of the four-particle interaction 
term by choosing the non-local inverse propagator 

P,(0 = ^[A7 + 2«:F,(/)]-\ (50) 

where A^^ is the microscopic (bare) contact coupling in Eq. ([8]). The modified bare action 
S'e now reads 

SE[^.r,x] = I rm^^ + f]^{Q) + \ [ x{-Q)PxiQ)x{Q) 

^9 I r{Qi)HQ2)x{Q3)S{-Qi + Q2 + Qs) 

JQi,Q2,Q3 

with the four-boson interaction replaced by the Yukawa-like term X"*/""^*- Our truncation for 
the flowing action is chosen to be 

T,[^l;,r,x] = [ r{Q)PAQMQ) + l [ x{-Q)PxiQ)x{Q) 

JQ ^ JQ 

+ g [ r{QiMQ2)x{Q3)S{-Qi + Q2 + Qs) (52) 

JQt,Q2,Q3 

r{QiMQ2)r{Q3MQ4)S{-Qi + Q2-Q3 + Q4), 



where A is a contact coupling, which is zero in the ultraviolet by construction and is regen- 
erated during the RG flow through a box diagram. 

Due to the numerous simplifications of the nonrelativistic vacuum, which are described 



in detail in 



33j . we note 
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the inverse propagators and are not renormalized during the RG flow 

2 

P4Q) = tco + q, P^iQ) = ^ [Aj:^ + 2kF,{1)]-' . (53) 

Technically this arises from the fact that the one-loop Feynman diagrams, which would 
renormalize P^ and P^, have poles in the same half plane of the complex loop frequency. 
One can close the integration contour such that it does not enclose any frequency 
poles. As the result of the residue theorem, the frequency integrals vanish, leading 
to the non-renormalization of the propagators. We observe the clear manifestation 
of the non-renormalization of the long-range potential, mentioned already in Sec. HTl 
Another example for this behavior is the boundary sine-Gordon theory, which was 



studied in 
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42|. 



The coupling g can be chosen arbitrarily, e.g. g = 1, at the UV scale. This is a well- 
know ambiguity of the Hubbard-Stratonovich decoupling. Furthermore, the coupling 
g is also not renormalized during the renormalization group evolution for the reason 
mentioned in the previous point. 

Interaction terms of the form with n > 1 and %l)*il)x^ with k > 1 are absent in the 
ultraviolet and are not generated during the flow due to the argument provided above. 
On the other hand interaction terms of the form {ip*ip)'^x^ with k > are generated 
through box diagrams during the renormalization group evolution. 

Our truncation is complete up to the two-body sector. This is a result of the special 



hierarchy of the flow equations in the nonrelativistic vacuum [32|, |33[ , where the cou- 
plings from the higher-body sectors do not influence the couplings from the lower-body 
sectors. 

It is straightforward to derive the flow equation of the momentum-independent coupling A. 
The result turns out to be given by Eqs. (IHl [15]) with the substitution F2{1) — )■ F^il) + 
in Eq. (fT5|) . The initial condition \^{k = A) = A^^ of the one-channel model of Sec. HT] is 
implemented directly as the initial condition of the inverse propagator P^ in the bosonization 
approach. Otherwise, the flows are completely equivalent at the level of our approximation. 
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B. s-channel bosonization 



We will follow an alternative bosonization procedure in this subsection. In this approach 
the four-particle contact vertex is approximated by the exchange of a massive (M^ = 2M^) 
particle in the s-channel. 

At the level of the microscopic (bare) action this is achieved via the Hubbard- Stratonovich 
transformation. The auxiliary part, added to the microscopic action ([8]) at T = /i = 0, reads 

SE,au.[^A\4>A*]= f rn\<p*{Q, + Q,) + -^r{Qi)r{Q3)}x 

{ct>{Q2 + Qa) + ^V^(Q2)^(Q4)}X (54) 

S{-Qi + Q2-Q3 + Q4). 

For the specific choice ^ = we can cancel the contact term in the modified microscopic 
action S'^ = Se + SE,aux 

JQ 

W{Qi)m2)m^) + (t>{QiW{Q2W{Q^)] K-Qi + Q2 + Q3] 



1,Q2,Q3 



[ Fa{i)r{QiMQ2)riQ3mQM-Qi + Q2-Q3 + Q4). 

(55) 

During the renormalization group evolution the massive field becomes dynamical and 
develops an inverse propagator P^{Q) = f{iuj + \ + m^), where / is some yet unknown 
function. The argument of the inverse propagator is fixed by the Galilean symmetry of the 
nonrelativistic vacuum, while the function / is determined by the dynamics. Our truncation 
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dth = at~ \ f + d, 




dt\ = dt 




+9, 




+dt 




+dt 



Figure 8: Flow equations of P(j){Q), h and A in the form of Feynman diagrams. The solid lines with 
one (two) arrow(s) denote regularized propagator for the field {<j)), the dashed lines correspond 
to the long-range interaction vertex, the dark blob represents the contact coupling A. 

of the flowing action within the s-channel bosonization approach is 



+ 



h 



{(t>*{QimQ2mQz) + miW{Q2W{Q^)) k-Qi +Q2+ qs) 



1,Q2,Q3 



- K I Fd{i)r{QiMQ2)r{Q3MQA)Si-Qi + Q2-Q3 + Qa) 

}i,Q2,---,Qi 

[ r{QiMQ2)riQ3MQA)si-Qi + Q2-Q3 + Qa), 

(56) 

where h denotes the momentum-independent Yukawa-hke couphng introduced above and A 
stands for the contact four-boson IPI vertex, which is regenerated during the flow through 
the box diagram. The non-trivial flow equations for the couplings of the average action f l56|) 
are depicted in Fig. [8] in terms of Feynman diagrams. 

It is desirable and possible to cancel the RG flow of the coupling A by absorbin g it into 



the flows of the other couplings. This can be achieved by a rebosonization procedure 
The idea is to make the composite fleld scale-dependent — 0^ with the choice 

where ak is some real, scale-dependent function. The one-loop flow equation 
average action generalizes to 
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44|. 



OkTk^-miTf^^Rk) 



kVk 



•k9k- 



(57) 



for the 



(58) 
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The unknown scale-dependent coefficient ak can now be fixed by the condition that the flow 
of the contact coupling A is zero 

= ^ A = 0. (59) 

The resulting flow equations (with unspecifled cut-off) for the Yukawa coupling h and the 
inverse propagator P(i,{Q) now read 

dtP^iQ) = -^/iVo,,(Q) 

. , An^P^iQ = 0) Jy 

dth = 2hKJiA H — ; 

h 

where Jn^ was deflned in Eq. f llSp . and for Jq^^ we introduce the momentum- dependent 
generalization 

Jod{Q)= — 7 — r^^t • 61 

J 2tt{2ttY [^^uj + uj,) + {l + q)^ + Rk{L + Q)][-zu + l'' + Rk{-L)] 

In order to proceed we must specify the cut-off function and the ansatz for the inverse 

propagator of the fleld 0. We choose the sharp cut-off, introduced in Sec. [TTl and take the 

sirnple ansatz P^{Q) = A^{iui + Y) +"^^ as used also for the related vacuum 3-body problem 



m 



34| . With this choice the flow equations become 



d,ml = (2{ - - 2)ml + 

a,.,.K-2--^).«-^^, (62) 

dtA^R = (2e - d)A^R - ^hl, 

where we expressed the flows in terms of the rescaled parameters = /c~^~'^+^^m^, A^ji = 
k~'^~^'^^A^, and hf> = The parameter ^ is not specifled yet, but can be 

identified with the naive scaling dimension of the composite field, i.e. [0] = ^. At this point 
it is convenient to multiply the second equation in ( l62l) by 2hR and obtain fiow equations 
for and hj^ 

d,ml = {2i-d- 2)m5; + ]-hl, 

The equations in f lB^ form a closed linear system of first order ordinary differential equations. 
The eigenvalues A± of this system can be readily found 

A± = 2e-3-^-^±^v/(rf-2)2-4«:. (64) 
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We note that = ^"^—p-, found in Sec. [TTl provides an important boundary value also in 
the current analysis. While for k, < n^r (repulsion and weak attraction) the eigenvalues are 
purely real, for n > Her (strong attraction) acquires a non-trivial imaginary part. 

Let us investigate the latter case in more detail. It is convenient to cancel the real part 
of by choosing = | (3 + | + ^^) . With this choice the eigenvalues A^*" are imaginary, 
leading to a pure oscillatory behavior of the flows of m|j and hj^. Choosing the initial 
conditions m'j^{t = 0) = M, hj^it = 0) = H, we obtain 



ml{t) = Mcos(^--5t) + W^Msin(^--5t) 

where D = {d — 2)^ — 4k, is the discriminant of the quadratic /3-function introduced in Sec. 
inland Wm and Wh can be expressed as 

2k 



W. 



M 



2-d 



d-2 



M + H 
2-d + 



d 



2k ^ ^''^ 
2k ^ 



\id-2y- 

As expected, we obtain a finite (bounded) limit cycle solution (!65|) for the pair (m|,, /i^) for 
strong attraction k > K^r- It is clear that the infinite number of the two-body bound states, 
present in this case, manifest themselves as zeros of the mass coupling Tn'j^{t) during the RG 
flow. 

It is possible to complexify the flow equations fl63|) in a similar manner as presented in 
Sec. IIVI by extending the couplings hj^ and m|. to the complex plane 

hi = hm + ihl^ m\ = + im\^. (67) 

In this extension the non-zero value of makes the composite particle (f) unstable. If 
small compared to m^^, it equals to the decay width, which is inverse proportional to the 
life-time of the particle. Due to the linearity of the evolution equations fl63|) . the flow 
equations for the pairs {fnli, h^^), {mj^^j^m) decouple. They have exactly the same form 
as Eq. flB51) with the solution flBSj ESD, provided we perform the substitution m|, — m^-, 
hf^ — hj^-jM — )■ Mi,H — i- Hi with i = 1,2. We depict the numerical solution of the 
complexified system in Fig. [91 As expected, the flow of the coupling A/? = 2^ reproduces 
our finding from Fig. [5l 
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Figure 9: Periodic flows of the real part (blue) and imaginary part (red) of the complex couplings 
^R' ""^R ^^"^ -^R ~ 2n?" specific choice of parameters d = 3 and k = 1. The amplitude of 

the m|j2 oscillations is related to the decay rate, which must be fixed by experiment. The width 
of the peaks of Xr2 is of the order of magnitude of the decay rate. 

Most importantly, we are now in the position to find an analytical solution of the general 
non-linear system of flow equations introduced in Sec. IIVI 



(68) 



with q;,/3,7 G M and D = jS"^ — ia'j < 0. In order to achieve this goal, it is sufficient to 
project the complex coupling 



"'R _ "-Rl ^ ''"'R2 



2ml 2[m|ji+im|2] 



(69) 



onto its real and imaginary part A/ji and Xr2 

_ ''^Rl^Rl + "^R2^/?2 \ 

— — TTi — 4 — ; — n — 



2[m|ji + m^2] '^[^Ri + ^m] 

This expression comprises an explicit solution of the system (!68|) provided we substitute 



(70) 



= Mi cos( " t) + WMi sin( t) 



(71) 



with D = (3'^ — 4a7 and 



{l3Mi + aHi 



\/-D 



(72) 



where z = 1, 2. 
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We comment on the issue of initial conditions. Obviously, the first-order system (168|) 
must be supplemented with two initial conditions Xniit = 0) and XR2{t = 0). In our 
description, however, there are four unknown constants Mi, M2, Hi and H2 to be fixed. 
The identification 

^^^^^ - °^ - 2 Mf + M| ^^^^^ - °^ - 2 Mf + M| 

fixes only two of them, leaving the remaining two ambiguous. This is the known arbitrariness 
arising in the Hubbard-Stratonovich transformation. Remarkably, one can fix the remaining 
two constants almost arbitrarily®. An especially convenient choice corresponds to Mi = 1 
and M2 = 0. With this choice we have 

Xm{t = 0) = ^Hi Xj,^{t = 0) = ^H2. (74) 

We present also the analytic solution of Eq. ( 168|) in the conformal phase with D > 0. The 
solution is still given by Eq. f l70|) provided we substitute 



with 



m\- = Mi cosh(— -t) + WMi sinh(— -t) 
hi, = Hi cosh(^t) + Wm sinh(^t) 



WMi = -4^{f3Mi + aHi) 

(76) 

Wm = -j= (47M, + l3Hi) . 




Finally we note that the analytic solution ( 169|) can be expressed in the unified, compact 
form 

e 2 + Ge 2 y 
for both conformal and non-conformal phase. In the last expression 

for D < 0, and the complex constant C determines the initial condition for Xr. Using Eq. 
(177]) we compute the curvature k and the related radius R of the complex RG trajectories, 
which are given by 



K 



(78) 



Some discrete values of two constants lead to a degenerate solution and are therefore forbidden. 
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Here the primes denote the first and second derivative with respect to t. We observe that 
the curvature does not depend on the "RG time" t in both the conformal and non-conformal 
phase. For D > the trajectories form arcs of circles of radius R = 1 2aiLc I (^^^ -^^S- E]), 
while for D < the trajectories constitute closed circles of radius R = \ a{i^\cp) \ V—D (see 
Fig. ED. 

VI. CONCLUSIONS 

In this work we investigated the nonrelativistic, quantum-mechanical problem of an in- 
verse square potential problem using functional renormalization. The potential is classically 
scale- invariant. Additionally, it is singular at the origin and the model must be augmented 
by a contact term, which is necessary for renormalization. We demonstrated that the RG 
flow of the contact coupling A either approaches a real fixed point (conformal phase) in the 
IR and UV or undergoes a scale anomaly which manifests itself in an infinite, unbounded 
limit cycle (non-conformal case). The overall behavior is determined by the sign of the dis- 
criminant D of the quadratic /3-function of the contact coupling, and depends on the strength 
of the long-range inverse square potential and spatial dimension. Remarkably, in the non- 
conformal phase {D < 0) the /3-function possesses a pair of complex conjugate fixed points. 
This observation led us naturally to the extension of the RG analysis to complex values of 
the coupling A — )■ Ai + 1X2- The RG evolution was computed numerically and the phase 
portraits of the flows were obtained. Additionally, we provided a geometric description of 
the complex flows on the Riemann sphere. We observed that in the non-conformal phase the 
real part Ai and the imaginary part A2 develop a finite (bounded) limit cycle in the complex 
plane. This suggests that the complex extension can be utilized as an efficient numerical tool 
for investigation of various problems with infinite (unbounded) limit cycles, as for example, 
in few-body physics of cold atoms near unitarity. One should simply add a small imaginary 
part A2 to the real part Ai and follow the extended RG evolution. This procedure regulates 
the periodic singularities during the RG flow making the numerical treatment feasible. From 
a physical point of view, the introduction of the complex coupling constant appears quite 
natural since it arises from an inelastic channel in the two-body scattering. The bosonization 
procedure allowed us to view the problem from a different perspective. More importantly, 
it enabled us to find an analytical solution of the extended set of non-linear flow equations. 
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Quite surprisingly, a similar RG behavior as in our case of an inverse square potential 



appears also in the context of chiral dynamics in QCD with a large number of flavors |45|, 146 1. 
Here, gluon exchange between the quarks induces quark self-interactions. The flow equations 
governing the RG evolution of these self-interactions exhibit the same behavior as shown in 
Fig. [2J Depending on the value of the gauge coupling a the flow equations possess either 
two fixed points (a < Ucr) and the system is in the chiral symmetric phase, or as a increases 
above this critical value acr one ends up in the chiral symmetry broken phase. The analog 
to the conformal symmetry breaking in dependence on the value of k in our simple case of 
an inverse square potential appears to be quite remarkable. 

Finally, we comment on a possible connection between our work and the recent studies 
of non- Fermi liquids using the AdS/CFT correspondence 47-49]. The authors of Q, 4^ 
studied the relativistic many-body physics of fermions at vanishing temperature in d space- 
time dimensions by mapping onto a classical gravity problem with an extremal charged 
black hole in anti-de Sitter spacetime {AdSd+ij- In this description the low-energy scaling 
behavior around the Fermi surface is related to the near-horizon geometry, which turns out 
to be AdS2 x R"'^"'^. Notably, the isometry group of the AdS2 part is 5*0(2,1), which is 
exactly the symm etry group of the quantum mechanics of the inverse square potential in the 
conformal phase 18|. This suggests that the emergent IR CFT, defined in 49|], might be 
conformal quantum mechanics with the inverse square potential. Another evidence in this 



the real and imaginary 
agrees remarkably 



direction is given by the observation that the low-energy behavior o " 
parts of the retarded Green functions, computed numerically in [48, 
well with the RG flows of the real and imaginary part of the contact coupling in conformal 
and non-conformal phases (see Figs. O [6]). This suggests that the complex extension and 
its analytical solution, found in this work, might be useful for a better understanding of 
non-fermi liquids from the AdS/CFT correspondence. 
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Figure 10: RG flows on the complex plane and the corresponding map onto the Riemann sphere 
in the non-conformal phase. The pink dots represent the two complex fixed points. The solid lines 
with arrows, pointing towards the UV, correspond to different RG trajectories. Additionally, the 
black straight solid line from the north pole illustrates the stereographic projection of the Riemann 



Appendix A: Complex RG flows on the Riemann sphere 

As has already been noted in Sec. IIVI the generalized RG flow of the complex coupling 
A = Xi + iX2 can be most conveniently studied on the Riemann sphere. The Riemann sphere 
is a unit sphere S"^ in three dimensional real space M^, which intersects the complex plane 
C = C U {00} at the equator. The map a: 5^ — t- C onto the complex plane (Ai,A2) is 
given by the stereographic projection (see Fig. \1U^ and reads in the Cartesian coordinates 
{x,y,z) G 



sphere S*^ onto the complex plane C. 



Ai 



X 




1-z 



(Al) 



The inverse map is given by 
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Figure 11: RG flows on the complex plane and the corresponding map onto the Riemann sphere 
in the conformal phase. The pink dots represent the two real fixed points. The solid lines with 
arrows, pointing towards the UV, correspond to different RG trajectories. 

2Ai 2X2 + A2 — 1 , , 

^= Af + Ai + 1' ^= A? + Ai + 1' ^= Af + Ai + 1' 

x = 0, y = 0, z = l, A = 00. (A2) 

The analytical solution of the complex flow equations (13T1) was obtained in Sec. |Vl We 
substitute these solutions into Eq. ( ]A2I) and obtain the analytical solution on the Riemann 
sphere, which we plot for both the non-conformal {D < 0) and the conformal {D > 0) phases 
in Figs. Uni and [IH 
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